Stabilizing Textures with Magnetic Fields. 



R S Ward* 

Department of Mathematical Sciences, 
University of Durham, 
Durham DHl 3LE 

February 1, 2008 



Abstract 

The best-known way of stabilizing textures is by Skyrme-like terms, but another 
possibihty is to use gauge fields. The semilocal vortex may be viewed as an example 
of this, in two spatial dimensions. In three dimensions, however, the idea (in its 
simplest form) does not work — the link between the gauge field and the scalar field 
is not strong enough to prevent the texture from collapsing. Modifying the |Z)<I>p 
term in the Lagrangian (essentially by changing the metric on the $-space) can 
strengthen this link, and lead to stability. Furthermore, there is a limit in which the 
gauge field is entirely determined in terms of the scalar field, and the system reduces 
to a pure Skyrme-like one. This is described for gauge group U(l), in dimensions 
two and three. The non-abelian version is discussed briefly, but as yet no examples 
of texture stabilization are known in this case. 
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1 Introduction 



Textures are classical solutions which are characterized by a nonzero homotopy group 
7rrf(T), d being the number of space dimensions. The relevant systems typically 
involve a scalar field $ taking values in the target space T. With a Lagrangian 
such as |9^^>p, and for d > 2, configurations are prone to implode (by the usual 
Derrick scaling argument). In an expanding universe, textures might be stabilized 
by the cosmological expansion; but we are interested here in cases where gravitational 
effects are negligible, and we take space-time to be flat. In flat space, the best-known 
way of stabilizing textures is to add a Skyrme term involving higher powers of dfj^^. 

By contrast, vortices or monopoles correspond to a nontrivial ird-iiT), and (in 
their "local" versions) are stabilized by gauge fields. Many similarities between tex- 
tures and vortices/monopoles have been noted. For example, multi-Skyrmions and 
BPS multi-monopoles (located at a single point in space) each have a polyhedral 
structure corresponding to an appropriate subgroup of 0(d), and this has been partly 
understood in terms of rational maps from the Riemann sphere to itself [Q] . The pur- 
pose of this paper is to investigate the stabilization of textures by gauge fields, and 
so in particular it explores a different sort of relation between the two classes of 
topological solitons, generalizing the example provided by the semilocal vortex [||, ^. 

The idea of stabilizing textures with gauge fields has been investigated before. 
One motivation has been the fact that Skyrme terms are non-renormalizable, whereas 
gauge theories may have better quantum behaviour; but in this paper the consider- 
ations are entirely classical. For the extended abelian-Higgs model (with the Higgs 
field being a complex doublet), it was pointed out in [Q] that an expansion in field 
gradients produced a Skyrme-like term, which suggested stability; at the time, this 
was not investigated in detail. More recent numerical simulations seemed to show 
that stability was indeed present (although, as reported below, we have not been able 
to confirm this result). In a different abelian system (involving a triplet of real scalar 
fields and a massive abelian gauge field) no stable textures could be found Q. For 
the non-abelian case, scaling arguments again suggest stability (c/. 0); but detailed 
investigation such as |p have produced negative results. The conclusion seems to be 
that the scalar field and the gauge field have to be linked to each other sufficiently 
strongly in order to prevent each from collapsing independently; and in "standard" 
systems, this link is not strong enough. 

The general framework is as follows. Suppose we have a system involving a gauge 
field (with gauge group G), and a multiplet <I> of scalar fields coupled to it. The 
"basic" Lagrangian of the system has the form 

C = i\D^<S>\^-i{F,,f-Vi<^). (I) 
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For space dimension d = 2, the system defined by (|l|) may admit stable static so- 
lutions (for example, semilocal vortices); but for d = 3 it seems not to — some 
modification is needed. The idea pursued here is that the term |D^<I>p in the La- 
grangian involves a choice of metric on the space T in which <I> takes its values, and 
we can change this metric. For example, if $ is a complex vector, then the standard 
Euclidean metric is iD^p = {D^^){D^), where D^^ denotes the complex-conjugate 
transpose of D^. A natural modification of this (see the following section) is to add 
a term K'^|<I>^Z)^<I>p, where k is a constant. So we now have a family of systems, 
parametrized by k. Taking the limit k ^ oo enforces the constraint 

^^D^<^ = 0, (2) 

which (under favourable circumstances) determines the gauge potential in terms of 
So we have a family of systems where, in an appropriate limit, the gauge degrees of 
freedom disappear, and the Maxwell/Yang-Mills term {F^y)"^ becomes a Skyrme term. 
This enables us to track a soliton solution as it changes from a gauge-stabilized texture 
into a Skyrme-stabilized texture. In an appropriate limit of parameters k, . . . — > oo, 
one gets a Skyrme system which certainly admits stable solitons; one question is for 
which finite values of these parameters there are stable solitons. 

Non-trivial examples of this idea have only been found in the abelian case G = 
U(l), and these are described in the next two sections (for d = 2 and d = 3 respec- 
tively). A discussion of the non-abelian case {G = SU(2) in d = 3) is given in the 
final section. The conclusion therefore is that textures can be stabilized by (abelian) 
magnetic fields, but no non-abelian version of this appears to be known. 

It might be noted that the idea of adding a term and investigating 

how solitons depend on the parameter k, has been investigated before; the simplest 
example (in a somewhat different context) is that of the CP^ model with no gauge 
field |]. 

2 Semilocal Vortices and Planar Skyrmions. 

In this section, we take d = 2 (so space is the plane R^), and gauge group U(l). 
Let the Higgs field $ be a complex doublet $ = [<I>^ <I>^]*. The resulting extended 
abelian-Higgs system admits semilocal vortex solutions [§, |3|; and in the limit k — > oo 
it becomes, as we shall see, a Skyrme version of the CP^ model. The generalization 
with $ being an M-tuplet, leading in the limit to a Skyrme version of the CP^^~^ 
model, is straightforward; but for simplicity we shall restrict here to the CP^ case. 
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The standard Lagrangian is 

£ = i (D^^)^ {D'^<!>) - i - iA(l - $t$)2 , (3) 

where -D^$ = 5^$ — iA^*^. For the semilocal vortex solution, the gauge field provides 
a "hard core" which prevents the soliton from shrinking. If < A < 1, the single 



soliton is stable; but for A > 1 it is unstable (it expands without limit) |10|, For 



A = 1 there is a one-parameter family of static solutions saturating a Bogomolny 



bound, but these solitons are marginally unstable |12]. One member of this family is 
(an embedding of) the standard Nielsen-Olesen vortex. 

Various relations between this system and the CP^ model have been noted before 
(c/. IQ, For example, imposing the constraint = 1 (this corresponds to 

letting the parameter A tend to infinity), and scaling away the {F^y)"^ term, leaves 



the CP^ model |13]. But in order to have stable semilocal vortices which become 
Skyrmions as a limiting case, one needs to make some modifications. 

Recall, first, the symmetry of this system The ungauged system has an 
S0(4) global symmetry. On gauging a U(l) subgroup, this S0(4) is reduced to the 
product of the local U(l) and a global SU(2); the field ^> belongs to the fundamental 
representation of this SU(2). The most general SU(2)-invariant metric on T = is 
hpqD^^ D^'^, where 

hpgi^,^) = g{C)5pQ+g{0^P^Q (4) 

and ^ = ^^^^ = ^tcj). The two functions g and g are arbitrary. But in the limit 
A ^ oo, which is of particular interest here, we have ^ = 1; so let us take g and g 
to be constants, scaling g to unity and writing g = . Using this modified metric 
instead of the standard Euclidean one amounts to replacing (|3|) by 

The second modification is as follows. In order to have stability for A > 1, we 
need an extra potential term, which necessarily breaks the SU(2) global symmetry 
(see for example [U, |l5|). We shall add to (^) the term a|$^p, where a is a positive 
constant. In the Bogomolny case (k = and A = 1), there is now a unique minimium: 
it has <I>^ = 0, and is the Nielsen-Olesen vortex with energy E = it. 

With these two modifications, the static energy density £ of the system is given 

by 

28 = {Dj^)^ {Dj^) + K^\^^Dj^^ + {Bjf + V{^) , (6) 

where V{^) = — $t$)2 _|_ 2a|$2p, and where Bj = ejkidkAi is the magnetic 
field strength. 
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The boundary conditions are chosen to ensure finite energy. At spatial infinity, 
one must have 

(a) A, = f-^djf, where |/| = 1; 

(b) Dj^ = ^ ^ = f^^K, where K is a constant 2- vector; 

(c) = 0=^ K =[k 0]* with \k\ = 1. 

Because of (b) and (c), ^ cannot be zero at spatial infinity; and in order for $ to be 
single- valued, / has to be single- valued. Hence / is a map from the circle at spatial 
infinity to the gauge group U(l), and the degree of / is the soliton number A^. The 
total magnetic flux is proportional to N, in the usual way. The fact that there is 
nontrivial topology does not necessarily mean that there are stable solitons; but the 
numerical work described below indicates that there are, at least for certain ranges 
of the parameters a, A and k. 

Taking the limit A ^ oo enforces the constraint <I>^<I> = 1 (so ^ takes values in 
S^). If in addition «; — > oo, then the minimum-energy configuration approaches one 
for which ^^Dj^ = 0, and hence 

Aj = -i^^dj^. (7) 

With Aj is given in terms of ^> by this expression, (Dj^)^ (Dj^) becomes the standard 
CP^ energy, and (Bj)'^ becomes a Skyrme term. We can re-express this as an 0(3) 
sigma-model in the usual way: define a unit 3-vector field ^ by i?^ = $^(T$, where 
a"' are the Pauli matrices. This corresponds to the standard Hopf map from (the 
space = 1) to S^. Strictly speaking, the $ field is a vortex (winding at spatial 
infinity); but the (j) field, obtained from it by projection, is a texture (constant at 
spatial infinity). Then in the A, k ^ oo limit, the energy density £ is given by 

8£x,^^oo = {dj$) ■ {dj$) + di$ X d2$f + 4a(l - n • 0), (8) 

where n is a constant unit vector. This is a planar Skyrme system pB, 17]. The 



energy of the Skyrmion solutions depends on a, and can be found by numerical 
minimization; for the 1-soliton with a = 1 it is = 3.15577r. 

The energy (and the stability) of the solitons in the system (|6|) may be investigated 
numerically, as a function of the three parameters a, A and k, and of the soliton 
number N . This has been done for the = 1 case, with a = 1 and A = \ + . 
The result is summarized in Figure 1, which shows the energy E as a, function of 
K > 0. It was obtained by assuming the standard form for 0(2)-symmetric fields, 
namely <I>^ = f{r)ex.p{iN9), $^ = g{r), Ar = and Ag = a{r), where /, g and 
a are real- valued. The discrete version of the energy functional E[f,g,a] was then 
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minimized numerically, using a conjugate- gradient method. For each value of k, a 
stable solution was found. Note that, as expected, E goes from E = tt (the Nielsen- 
Olesen vortex) at k = and A = 1, to ii^ = 3.15577r (the planar Skyrmion) as k ^ oo 
and A — > cxd. 
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Figure 1: The energy E of the 1-soliton on R^, as a function of k, with a = 1 and A = 
The dashed line is the energy of the planar Skyrmion obtained in the limit k oo. 



3 Vortex Rings and Hopf Textures. 

In this section, we investigate the same system (^ as before, but in spatial dimension 
d = 3. The extra potential term is omitted (in other words, a = 0), so the global 
SU(2) symmetry is unbroken. One may form a texture configuration by taking a 
finite length of semilocal vortex with its ends joined together to form a loop in 3- 
space, and it has previously been speculated that such a texture might be stable 

ii- 

The energy density is given by @, with a = 0; so the system depends on the 
two parameters A and k. In the limit A,k ^ oo, we again get an S'^-valued scalar 
field $, with the gauge potential being given by (^); it has previously been pointed 
out (c/. [|^, |l^) that this limit is equivalent to the Faddeev-Hopf system 22, 
|23| , p^ , p5| , [2^ ]. So there are stable ring- like solitons in the limit; the question here is 
whether they are stable for finite values of k and A. 

The boundary conditions imply, as before, that <I> = [<I>^ <I>^]* = f^^K at r = 00, 
where K is a, constant 2- vector; so W = /^'^ is constant at spatial infinity. Think- 
ing of as a stereographic coordinate for CP^ therefore shows that $ (provided 
it is nowhere-zero) defines a map from R'^ to S"^ which is constant at infinity, and 
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hence is classified topologically by the Hopf number G 713(6'^). For = 1, the 
field resembles a single vortex ring. The stability of such = 1 configurations 

1.74 1 , , , , , , 1 
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Figure 2: The energy E of the 1-soliton on R^, as a function of k, with X = 1 + k'^. 



has been investigated numerically, again by minimization of the energy functional. 
The solitons cannot be spherically-symmetric, but one expects that for small values 



of N they will be axially-symmetric |21, 22|. So one can reduce the problem to a 
two-dimensional one which is not too difficult computationally. More precisely, one 
can use cylindrical coordinates, and impose an S0(2) xS0(2)-invariant ansatz, as for 
example in [^. 

Minima were sought for the one-parameter family of systems obtained by setting 
A = + 1, and stable solitons were found for k > 7.1. Their energy is plotted in 
Figure 2. For n < 7, however, the radius of the vortex ring shrinks to zero, and 
the field unwinds: there is no stable minimum. When k (and therefore A) tend to 
infinity, the normalized energy E' = E/Att approaches the value E'^ = 1.73 (obtained 
by extrapolation of the data in Figure 2). This is exactly the energy of the single 
Hopf soliton: from and allowing for different coupling constants, we get the 
value E' = 1.22^2 = 1.73. 

The analogous computation previously reported in Q for the k = case suggested 
that one might have stability for fairly small values of A (of order unity) . The results 
described above do not confirm this, and in fact no stable solution could be found 
for K = 0, even with A quite large. But (as emphasized in |5|), there might be local 
minima in the configuration space which are difficult to detect, and which require 
an initial condition which is very close to the actual solution. So it remains an open 
question as to whether stable vortex rings exist for small values of k and A. It is, 
however, the case that the configuration which is stable for k = 7.1, A = 51.4 collapses 
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if K and A are reduced to k = 7, A = 50. 



4 Non-Abelian Gauge Field. 

As mentioned in the introduction, the question of whether textures can be stabiHzed 
by a non-Abehan gauge field has previously been investigated; there are suggestions 
based on simple scaling arguments (c/. Q), but more detailed studies have yielded 
negative results (c/. §). Let us look at the three-dimensional case {d = 3), with 
gauge group SO (3). The field $ belongs to some representation T of SO (3); so we 
have to choose F, as well as an appropriate potential function V{^). For example, 
for the 'tHooft-Polyakov monopole one uses the fundamental representation F = 3. 
The simplest extension of this is the four-dimensional representation F = 1 © 3. 
The corresponding system admits monopole-like soliton solutions which have been 
referred to as semilocal monopoles Q. (Another simple extension is F = 3 © 3, the 
corresponding solitons being referred to as coloured monopoles 

Let us look at the 1 © 3 case: so <1> = {(j)o,(j)) is a four- vector. Take the po- 
tential function to be A(l — so for large A, we get the constraint |<I>p ~ 1. 
One may then impose texture boundary conditions (rather than monopole boundary 
conditions): namely, $ tends to a constant as r — > oo in R^. So ^ is effectively a 
map from to S^, and it has a winding number N. The stability of spherically- 
symmetric = 1 configurations has been studied numerically — the details are as 
follows. 

For simplicity, we shall take the A — > oo limit, so = 1; and the metric on 
$-space to be flat (no extra term analogous to «;^|$^Z)j<^p). The energy density is 

£ = l\D,^f + ^,{Fjk?, (9) 

where Dj^ = {djcjiQ^djcjf — 2e°''"^A^(j)'^). To implement spherical symmetry, we take 
$ and the gauge potential A to have the standard "hedgehog" form 

A'^ = ejakx''f{r)/r^, ^o = cosg{r), = x'^ sin 5(r)/r , (10) 

with the boundary conditions /(O) = 0, /(oo) = |, g[0) = it, g{oo) = 0. The energy 
density then becomes 

^ = | + ^^^4 + ^[i + W-i)l. (U) 

One can then minimize the energy numerically; this was done using a conjugate- 
gradient method, with various initial conditions. But no smooth minimum could be 
found — in every case, both / and g collapse towards being zero almost everywhere. 
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One can see this collapse analytically, in the following highly-simplified version 
(involving just two degrees of freedom a and Let a and /3 be the values of r such 
that /(/?) = 1/4 and = 7r/4. In other words, a and /? are the "radii" of the 
scalar field and the gauge field respectively. More explicitly, take / and g to have 
the form 



One can compute the energy S(a, /?) of this configuration exactly: it is a rational 
function of a and /?. In particular, for /? = 1/a the energy has the form E{a^ = 
a X (polynomial in a). The salient point about this form is that its minimum occurs 
when a = 0; this corresponds to the scalar field shrinking to zero width, while the 
gauge field spreads out. As one sees from the usual Derrick scaling argument used in 
0, the contribution to the energy from the term can be reduced by scaling 

one way, while the contribution from {Fjk)^ can be reduced by scaling the other way. 
But the system as a whole can never reach a balance: the gauge field and the scalar 
field are just not sufficiently strongly coupled to each other to prevent each one from 
collapsing separately. As was remarked before, these results do not actually prove 
the absence of a stable solution: there might still be a local minimum somewhere in 
the configuration space. But it seems rather unlikely that this system does admit a 
stable texture. 

As in the previous sections, one can modify the metric on $-space, and this may 
improve the stability properties. That possibility has not yet been fully investigated; 
but certainly there is no gauge-invariant extra term, the vanishing of which deter- 
mines the gauge potential as in the abelian case (0). So the idea of obtaining the 
usual Skyrme model as a limit does not quite work in this non-abelian case. Using 
other representations F, and for that matter other gauge groups, opens up many 
more possibilities, which are still to be explored. But for the time being, it remains 
the case that there are no known examples of three-dimensional systems in which a 
texture is stabilized by a non-abelian gauge field. 
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